Transport on weighted Networks: when the correlations are independent of the degree 
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Most real-world networks are weighted graphs with the weight of the edges reflecting the relative 
importance of the connections. In this work, we study non degree dependent correlations between 
edge weights, generalizing thus the correlations beyond the degree dependent case. We propose a 
simple method to introduce weight- weight correlations in topologically uncorrelated graphs. This 
allows us to test different measures to discriminate between the different correlation types and to 
quantify their intensity. We also discuss here the effect of weight correlations on the transport prop- 
erties of the networks, showing that positive correlations dramatically improve transport. Finally, 
we give two examples of real-world networks (social and transport graphs) in which weight-weight 
correlations are present. 
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I. INTRODUCTION 

Complex networks have proved to be useful tools to 
explore natural or man-made phenomena as diverse as 
the Internet human societies [1], transport patterns 
between airports^, 0| or even metabolic reactions in the 
interior of cells [5|. The vertices in the networks repre- 
sent the elements of the system and the edges the inter- 
actions between them. The study of the topology of the 
network provides valuable information on how the ba- 
sic components interact. While the existence or not of 
an edge is already informative, in many cases, as those 
listed above, the interactions can appear on different lev- 
els. The bandwidth between two servers on the Internet 
for instance is not a flat quantity equal for all pairs, it 
depends on the importance of the servers as well as on 
the traffic expected. This fact led to the introduction of 
weighted graphs as a more accurate way to describe real 
networks [6|, |7i] . Weighted graphs are complex networks 
where the edges have a magnitude associated, a weight. 
The weight accounts for the quality of a connection. The 
existence of a distribution of weights dramatically alters 
transport properties of networks like the geometry of the 
optimal paths [1, [l^, [ll|, the spreading of diseases 
[12j or the synchronabizability of oscillators _13]. Most 
previous studies have been carried out on networks with 
uncorrelated weights on neighboring edges (those arriv- 
ing at the same node) even though most real cases possess 
correlation. Our aim here is to check how the presence 
of correlations can influence these results. 

There may be several kinds of correlations in random 
graphs [3]. Recently, it has been shown that the edge 
weights in some real- world networks are related to other 
properties of the graph such as the degree (the number of 
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connections a vertex has) [T^, [l^, ■ The weights were 
found to follow, on average, a power law dependence on 
the degree. Several theoretical mechanisms have been 
proposed to generate networks of this type [3. In this 
clear correlation is introduced between the weight 
of neighboring edges but one may wonder whether this is 
the only possibility for weight correlations. If not, which 
other structures are possible? How can the correlations 
be quantitatively characterized? And most importantly, 
which influence do they have on the transport? 

In this work, we address these questions. The organi- 
zation of the manuscript is the following. In the Section 

II, we present a model that allows to explore the differ- 
ent configurations for weight correlations independently 
of other properties of the network. Next, in the Section 

III, we consider and evaluate different magnitudes to esti- 
mate the type and intensity of weight-weight correlations. 
Section IV includes a study on how the presence of weight 
correlations affects transport. In Section V two examples 
of real-world networks showing this type of correlations 
are discussed: the IMDB actor collaboration network and 
the traffic network between US airports. And finally, we 
conclude and summarize in Section VI. 



II. A SIMPLE METHOD TO INTRODUCE 
WEIGHT- WEIGHT CORRELATIONS 

Let us start by defining a mathematical framework for 
the weight correlations. From the point of view of an 
edge of weight w with vertices with degree k and k' at its 
extremes, the joint probability that its neighboring edges 
have a certain weight is given by 

Pkk'{w,Wi,. . .,Wk-i,w[, . . . , Wfe,_i). (1) 

These functions contain all the information about both 
degree and weight distributions and correlations. How- 
ever, a situation in which a full hierarchy of such func- 
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FIG. 1: (Color online) Two possible cases in networks with 
correlations in the link weight: a) positively correlated nets 
and b) anticorrelated networks. The width of the line of the 
links represents the value of the weight. 



tions were needed to characterize the network would be 
hard to control from an analytical or numerical point of 
view. Therefore we will focus here only on the simplest 
scenario. In the same way the Markovian condition is a 
simplifying assumption for stochastic processes, we will 
consider only correlations generated by two-point joint 
probability functions Pkk'(w, w'), and, among those, ini- 
tially only the ones that are degree independent given by 
functions of the type P{w,w'). 

In order to construct weighted networks along these 
lines, we use the so-called Barabasi- Albert (BA) model 
[l9| . where new nodes entering the network connect to 
old ones with a probability proportional to their degree 
[20| . The networks generated by this model are scale-free 
(their degree distribution goes as Pk{k) ~ k~^), have 
no degree-degree correlations, and their clustering co- 
efhcient (probability of finding triangles) tends to zero 
when the system size tends to infinity. All this makes 
them ideal null models to test correlations between edge 
weights. Once the network is grown, a joint probability 
distribution for the link weights P{w, w') and an algo- 
rithm for weight assignation are needed. With the func- 
tion P(w,w') one can calculate the weight distribution 
P['w) = J dw' P{w,w'), and the conditional probability 
of having a weight w' provided that a neighboring link 
has a weight w, P{w'\w) = P{w,w')/ P{w). We start 
by choosing an edge at random and giving it a weight 
obtained from P(w). Then we move to the nodes at its 
extremes and assign weights to the neighboring links. To 
do this, we follow a recursive method: if the edge from 
which the node is accessed has a weight wq, the rest, 
wi, . . . , Wk-i, are obtained from the conditional distribu- 
tions P{wi\wi-i). The recursion is necessary to increase 
the variability in case of anticorrelation (see below). If 
any of the links, j, has already a weight, it remains with- 
out change and its value affects the subsequent edges 
j + 1, . . . — 1. We repeat this process until all the edges 
of the network have a weight assigned [2l[ . 

For P(w,w'), we have considered different possibilities 



but here we will focus only on the following three: 



Pu{w,w') = ^^j^^, (2) 

P-{W,W') = (^^,+-1)1+. , 

where X+ = 2"a(l + a), Xu = and X_ = 
c^^/2^i(q^j a, 1 + a, —1) are the normalization factors for 
the distributions on the domain of weights (1, oo), and 
2Fi() is Gauss hyper-geometric function [24|. With- 
out losing generality, we have chosen these particular 
functional forms due to their analytical and numeri- 
cal tractability. The distributions generated by Eqs. 
([2]) asymptotically decay as P{w) ^ The rea- 

son to use power-law decaying distributions is that em- 
pirical networks commonly show very wide weight dis- 
tributions that in a first approach can be modeled as 
power-laws (see Fig. 6 and Refs. [1, i, H, HI). We 
name the functions as -I- (positively correlated), — (an- 
ticorrelated) and U (uncorrelated) because the average 
weight, {w){wq) = J dw w P{w\wo), obtained with the 
conditional probabilities from a certain seed wq grows as 
(w)-\-{wo) — {l+a+wo)/a, decreases as {w)-{wo) = {ce+ 
l/wQ)/{a — 1) and remains constant {w)u = a/{a — 1), 
respectively. This means that in -I- networks the links of 
each node tend to be relatively uniform in the weights 
(see Fig. la), with separate areas of the graph concen- 
trating the strong or the weak links, while in the negative 
case links with high and low weights are heavily mixed. 

From a numerical point of view, we have checked how 
the variables to measure vary with the network size N. 
In the following, most results are shown for N — 10^, 
which is big enough to avoid significant finite size effects. 
For each value of the exponent a (from Eqs. ^) and 
for each type of correlations, we have averaged over more 
than 600 realizations. Note that we use a as a control 
parameter for the strength of the correlations. For high 
values of a, P{w) decays very fast and the correlations 
become negligible, all links have almost the same weight. 
When a decreases however, the higher moments of P{w) 
diverge and one would expect the correlations to be more 
prominent. 



III. MEASURES OF WEIGHT CORRELATIONS 

After a look at the sketch of Fig. 1, the first estimator 
to consider in order to estimate weight correlations is the 
standard deviation of the weights of the links arriving 
at each node. If the weights are relatively homogeneous, 
the standard deviation will be lower compared with its 
counterpart in a randomized instance of the graph. The 
opposite will happen if the correlations are negative as in 
the case b of Fig. 1. More specifically, for a generic node 
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FIG. 2: (Color online) The rates E^, (diamonds), Tq, and 
p (circles) for the "+" distribution model as a function of 
a. The highlighted area encloses the resolution limit of the 
estimators. The rate is shown for several values of g: q — 2 
(squares), q — 4 (up triangles) and q — 6 (down triangles). 



of the network i, <7io{i) can be defined as 

jGi/(-i) 



(3) 



where is the set of neighbors of i and is the 
mean value of the weight of the links arriving at i. 
Once the deviation is calculated for each node, an av- 
erage can be taken over the full network getting (cr^,) — 
(l/A^) cru,(i). Then to evaluate the effects of weight 
correlations, it is necessary to compare the value of 
{o'w)org obtained for the original network with that mea- 
sured on uncorrelated graphs. It is, of course, important 
that the statistical properties of such uncorrelated graphs 
are as close as possible to those of the original graph. The 
most accurate procedure consists in disordering only the 
weights of the links of the original network. To do so, we 
interchange the weight of each link with that of a ran- 
domly selected edge preserving so the weight distribution 
P{w) and the network topology; z.e., degree distribution, 
degree-degree correlation, clustering, etc, remain invari- 
ant. Once (ctu,) is estimated for the original graph and 
for an ensemble of weight-reshuffled instances of it, the 
rate 



(4) 



can be calculated. If S^, > 0, the weight correlations in 
the original graph will be as in the case b of Fig 1. If 
it is identically one, there will be no weight correlations 
and if E^, < the correlations will be as in Fig. la. The 
behavior of S^j for the positive and negative models pro- 
posed in the previous section is displayed in Fig. 2. The 
first thing to remark is that indeed E^, can distinguish 
between the three cases. Moreover, it provides a first 
method to quantitatively estimate of the intensity of the 
weight correlations. 



FIG. 3: (Color online) Plot of the disparity as a function of 
the degree, the triplets represent the graphs generated using 
the "+" (blue diamonds), the uncorrelated (red squares) and 
"-" (black circles) distributions. Each triplet correspond to 
a different value of a, from bottom to top a = 20, 1.5, 1.2 
and 0.5. The straight line has slope —1 and is meant as a 
guide to the eye. The triplets have been shifted upwards to 
facilitate visibility. The inset shows the disparity versus k 
for the U case and for the values of a from bottom to top 
a = 20, 2,1.5,1.2,1 and 0.5. 



A similar result can be obtained with a magnitude 
that was previously discussed in the literature [1, |2^ . 
Its name is disparity and was introduced in the context 
of weighted graphs by Barthelemy et al. as a way to 
estimate how homogeneous the weights of the links ar- 
riving at a vertex are. The generalized disparity of node 



Yq{i), is defined as 



(5) 



where Si is the strength of i, Si — J2jeu{i) "^ij- 
links of a node have a similar weight, their value will be 
w « s/k, and therefore the disparity decays as Y2{k) ~ 
l/fc"^"^. On the other hand, if the vertex strength is 
essentially due to the weight of a single link, Yq{k) will 
tend to a constant. Typically instead of a generalized Yq, 
most of the works in the literature has focused on Y2, for 
which is commonly reserve the name of disparity. This 
latter magnitude can be related to by the following 
expression for each node i of the network 



Y2{^) 



1 

hi 



(6) 



where {'w)i is the average weight of the links of i and ki 
its degree. An important question to mention here is that 
the profiles of Yq{k) depend on the weight distribution, 
even for completely uncorrelated graphs. Yq{i) provides 
information on how different the weights of the links of 
i are but not on whether that particular configuration 
is or not a product of randomness or correlations. The 
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FIG. 4: (Color online) In the main plot, the dependence of p 
on the exponent a for the positively correlated graphs gener- 
ated with the "+" distribution of Eq. ([2J. In the inset, the 
same function but for the graphs of the "-" distribution. 



variation of Y2{k) with the exponent of the weight dis- 
tribution a for the uncorrelated U -model can be seen in 
the inset of Figure 3. In the same Figure, we also show 
the behavior of Y2{k) for the other two models, and 
" — " , for a few values of the exponent a. As before, in or- 
der to estimate the importance of the weight correlations, 
the disparity of the correlated graph has to be compared 
with that obtained from uncorrelated networks. To ex- 
press this comparison in a single rate, we can take the 
average of the disparity over all nodes of the network, 
(Yq), and calculate 



{^q} orig 
{^q} rand 



(7) 



The value of for the positive model as a function of 
a is displayed in Fig. 2. This magnitude is also able to 
discriminate between the different correlations. 

However, both E and Tq have resolution problems. As 
can be seen in Figure 2 for the positive model, if an area 
enclosing the numeric error is set immediately below one, 
the estimators S, Tq fall in turn relatively fast in that 
zone. The weights of the links in the model are 
continuous variables and therefore they are always cor- 
related. Although, as explained before, for higher values 
of a the effects of weight correlations can be weaker but 
still until a is not infinite they are not zero. An ideal es- 
timator should be able to distinguish the -I- model from 
a complete random configuration at very high values of 
a. In this context, seems to be the worst estimator. 
Tq is better than E^ and improves the higher q becomes. 
The reason for this behavior is that these magnitudes are 
not only estimating how wide the spectrum of values of w 
for a node is, they also supply information on the shape 
of the distribution of those values. As an example, let us 
consider a node with k links. The value of cr^ is higher if 
/c — 1 of them have weight a and the remaining weight 6, 
ayj = \a — b\ \/k — Xjy/k^ than if the distribution is more 
symmetric, let us say, with half of them with w = a and 



the other half with weight b, = \a — b\/{2^/2). The 
goal here is to study how different the amplitude of the 
weight values is compared with a random configuration 
of weights, hence the extra information contained in 
or Yq can be neglected. An ideal estimator for weight 
correlations only needs to consider the interval \a — b\. 
Following this idea, we define the range for a node i as 



(8) 



where Wmaxii) and Wmin{i) are respectively the maxi- 
mum and minimum weights of the edges of i. The de- 
nominator is a normalization factor to keep between 
zero and one. Note that r has a similar behavior to Yq in 
the limit g — > oo: if all the weights are equal r — Q, and 
also Yq ^ fc^~' ^ if g ^ oo. On the other hand, if the 
weight of link dominates the others r — > 1, ^ 1 too if 
q ^ CO. As before, to generate a correlations estimator, 
the average of r^, (r), can be taken over all the nodes 
of the network and contrasted with the equivalent value 
obtained from a set of weight-reshuffled instances. We 
will call p to the rate between these two quantities, 



(9) 



If p < 1, the network displays positive weight correla- 
tions. The stronger they are, the smaller p becomes. 
Otherwise, if p > 1, the weights are anticorrelated. 
p=T = E = lis the limit of uncorrelated networks. 
As can be seen in Fig. 2, p is a much more acute esti- 
mator of weight correlations than or E. Hence, from 
now on we will present our results as function of it. The 
variation of p with a is displayed in Fig. 4 for the 
and " — " models. The intensity of the correlations for the 
" -|-" model grows when a — > (p+ decreases for smaller 
a), while for the negative case p_ initially grows, peaks 
around a « 2.5 and then tends to one again for smaller 
a. 

Another question that we are now in position to board 
is in which way a relation between weight and degree af- 
fects the weight-weight correlations. As mentioned in the 
introduction, networks of this kind may display quite dif- 
ferent tr ansp ort pro perties from their unweighted coun- 
terparts O m, in, HI] . Usually, the weights of these 
networks are obtained by means of a relation of the kind 
Wij ~ {ki kjY [l^. The result is that provided that the 
degree is an " a-priorr characteristic that equally influ- 
ence all the edges of a vertex, the weight of the links is 
positively correlated. The networks created in this way 
show correlations similar to our model (regardless 
of the sign of the exponent 6) [1^. For instance, for 
S = ±0.5 the value of p is in both cases p « 0.832. 



IV. TRANSPORT PROPERTIES 

Let us focus now on the transport and how it varies 
with the presence of weight-weight correlations. Several 
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FIG. 5: (Color online) Relation between the intensity of the 
correlations, estimated with p, and the change in the weight 
of the superhighways for the positively correlated networks 
generated with the "+" model of Eq. ([2]). The values of a go 
from Of = 1 for the first point on the left to a = 2.4 in middle 
changing regularly by 0.2 for each point. The last point of 
the right hand side is for a = 20. 



measures have been proposed to study transport [13, [13] ■ 
In this work we will concentrate on the size and weight 
of the " superhighways^^ as recently introduced by Wu et 
al. [l^ . If the edges of an uncorrelated network are sev- 
ered following an increasing order from small to higher 
weights, the percolation threshold is eventually reached. 
The remaining connected graph after the process is over, 
the incipient percolation cluster (IIC) or superhighways, 
holds most of the traffic of the original network. Regard- 
ing this procedure, there are several questions to men- 
tion. First, along this work we have assumed that higher 
weights imply better transport properties. This is also 
the case for the two empirical examples discussed below; 
the weights represent number of passengers in one net- 
work and number of collaborations done together in the 
other. However, for some other graphs, the weights can 
mean higher resistance to transport. And, therefore, the 
superhighways must be obtained with the same proce- 
dure but cutting the links with highest weights first. In 
such circumstances, the lower the total weight of the su- 
perhighway is, the better the transport results. 

Another question is that the percolation threshold de- 
pends on the topological characteristics of the network. 
For uncorrelated undirected graphs, it is attained when 
in the process of severing random links, the rate (k'^) /(k) 
of the remaining graph approaches two, (fc^) / (k) =2 [27| . 
For directed networks or networks with high clustering, 
the threshold does obey different expressions [1^ . In the 
next section, we will further discuss this point as well as 
the numerical determination of the percolation threshold 
for our two empirical networks. 

Finally, it is worthy noting that this method to esti- 
mate superhighways cannot be applied directly to net- 
works with weight correlations. If the high weight links 
are concentrated in some areas of the graph, cutting the 



weak links is not going to affect those areas. This can 
change the percolation criterion, i.e. {k"^) / (k) = 2 can- 
not be used, although the percolation threshold that de- 
pends solely on the topology remains the same as in an 
uncorrelated graph. 

The goal in our case is to compare graphs with and 
without weight correlations and to quantitatively esti- 
mate the effects of these correlations. The method used 
in practice is to disorder the weights of the links of each 
correlated network. Then we estimate the superhighways 
of the randomized graphs and measure how many edges 
must be cut on average to attain the percolation thresh- 
old. Reac hing the percolation may have some numerical 
problems [2^ so the process of reshuffling the weights 
must be repeated several times. Next, we cut the same 
number of links in the correlated network (again going 
from lower to higher values of w) and compare the size 
and the weight of the biggest remaining connected cluster 
(Wsphgw) with the average of those found for the random- 
ized graphs. In this way, we obtain 

^'■sphgw — (W,pfcg„(rand)) ' 

g _ Number nodes in superhighway (orig) 
sphgw (Number nodes in superhighway (rand)) ' 

(10) 

The results, displayed in Fig. 5, show that, in general, 
positive correlations play a decisive role on the value 
of W sphgw, increasing it by orders of magnitude. The 
smaller is p respect to the unit, the stronger the effect of 
the weight correlations on the transport becomes. This 
phenomenon may be understood by keeping the analogy 
with the roads: the transport improves if the highways 
are connected together forming a communication back- 
bone as large as possible. Anticorrelated networks, on 
the other hand, exhibit smaller superhighways than their 
randomized counterparts although the effect is subtle. 

So far we have discussed models for which the weight 
correlations are independent of other structural factors. 
In general, there may be other aspects influencing the 
transport properties of a graph. If several compete, as it 
happens in the case Wij = {ki kjY , with S < 0, between 
the degree and the weight of the connections of a node, 
the transport capability of the network may suffer. For 
example, we measure ilsphgw — 75(1) for networks of 
size N ^ 10^ and (5 = 1/2, while ^sphgw = 0.014(1) if 
<5 = -1/2. 



V. EMPIRICAL NETWORKS 

Finally, we also consider a couple of real-world exam- 
ples. First the IMDB movie database with 383640 actors 
that are connected together whenever they have shared 
a common movie [l^, [s^]- This network is formed by 
the union of cliques, which means that the number of 
links is high, 15 038 083 in total. The weight of each 
link represents the number of times a partnerships has 
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FIG. 6: (Color online) In a), cumulative weight distribution 
Cw = f dw P(w) for IMDB actor collaboration network. In 
the bottom, in b), the clustering as a function of the degree. 
The line below corresponds to the function l/(fc ~ 1). 



been repeated. Higher values of the weight imply an in- 
creased probability of information transfer between two 
individuals. The cumulative distribution of weights for 
this network can be seen in Figure 6a. It shows a very 
wide functional form that can be well represented by a 
power-law. The presence of weight correlations in collab- 
oration networks have been discussed using a different 
technique in Ref. Here we will focus on the results 

obtained with p, ^sphgw and Ssphgw 

First of all, it is important to note that collaboration 
networks typically do not show a relation between the 
weight of the links defined in this way, or as social close- 
ness 7], and the degree of the nodes [3l|. Hence the 
weight correlations, if exist, are not a product of this 
type of relation. And, indeed, they exist since the actor 
network presents a value of p = 0.268(1). 

The measure of the superhighways in this case poses 
a certain level of challenge. The actor collaboration net- 
work, as many social networks, presents high clustering. 
The clustering of a node i is defined as 



2s,: 



(11) 



kiiki 1) 

where Si is the number of connections between the neigh- 
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FIG. 7: (Color online) Relative size of the giant component 
as a function of the removal probability q for the actor col- 
laboration networks. In the inset, the same quantity but 
in log-log and versus g* — q. Three values of q, are dis- 
played: (J, — 0.9987 (triangles up), qc = 0.9990 (circles) and 
q, — 0.9993 (triangles down). 



bors of i and ki stands for its degree. The average over 
all the nodes of the network with the same degree k can 
be then taken to obtain c{k), which is depicted in Figure 
6b for this network. In the same plot, it is also included 
the curve l/(fc — 1). The comparison is necessary because 
it has been recently shown that the percolation threshold 
of a graph is highly dependent on the clustering [2^ . In 
fact, clustered networks can be classified into two major 
groups: those with weak and those with strong cluster- 
ing. The difference between the two groups is whether 
c{k) decays as l/{k — 1) (weak clustering) or in a slower 
way (strong). The actor network clearly falls into this 
latter group. 

For weak clustered networks, it is possible to find a gen- 
eralization of the percolation threshold condition men- 
tioned in the previous section {{k'^)/{k) = 2) [28j. How- 
ever, as far as we know, nothing similar has been pro- 
posed for strong clustered graphs. Therefore, we will be 
force to use a more pedestrian technique to estimate the 
percolation threshold of the actor collaboration network. 
In Figure 7, the behavior of the rate between the size 
of the giant component gcc and its original value gcco 
is displayed as a function of the percentage of links sev- 
ered q. A continuous transition can be observed with 
this rate as order parameter. Assuming functional form 
of the type gcc/gccQ ^ {qc — q)^ , we find that the crit- 
ical point happens at a remotion rate of qc w 0.9990(3) 
(see the Inset of Fig. 7). The point in which the con- 
dition {k'^)/{k) = 2 is fulfilled, for instance, lies in a 
smaller value q = 0.9976(1). Once qc has been mea- 
sured, we can proceed as in the previous section, cut- 
ting a fraction qc of links following an ordered sequence 
from lower to higher values of the weight and compar- 
ing the results with those obtained for a graph in which 
the weights of the links have been reshuffled. The rates 
for the total weight of the superhighways for a few val- 
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FIG. 8: (Color online) In the main plot, the cumulative weight 
distribution for the US airport traffic network. The average 
weight of the outgoing connections is displayed in the inset as 
a function of the out-degree. 



ues of q are ^sphgw{q = 0.9976) = 3.2(1), flsphgw{q = 
0.9987) = 17(5), and flsphgwiqc = 0.9990) = 106(30). 
As can be seen, finding the value of Qsphgw requires a 
fine determination of q^- Even so, the high values of 
this rate gives us a clear feeling of the importance that 
the weight correlations have on the transport properties 
of these real world graphs. The values of Sgphgw that 
we find for the same removal rates q are Ssphgw{q ~ 
0.9976) = 0.62(2), Ssphgw{q = 0.9987) = 3.2(1.5), and 
Ssphgwiqc = 0.9990) = 20(7), respectively. 

The second example is a network composed by 1278 US 
airports. A directed edge connects two airports whenever 
there is a direct flight between them. The weight of the 
links represent in this case the number of passengers on 
that traject during 2005 The cumulative weight dis- 
tribution of this network is displayed in Figure 8. The 
weight distribution is in this case also wide, although 
clearly it does not follow a power-law decay. As can be 
seen in the Inset of that Figure, the average weight of the 
outgoing links exhibits a dependence on the out-degree of 
the nodes, {wij) ~ {kikj)^ '^'^ . Therefore it is not strange 
that the value of p that we measure, p = 0.983(1), de- 
lates the presence of positive weight correlations. Since 
the number of passengers in each direction can be dif- 
ferent, to calculate the superhighways it is necessary to 
generalize the concept to directed graphs. This means to 



study the incipient strongly connected component (SCC) 
instead of the incipient percolation cluster. Applying the 
same technique as the one illustrated in Fig. 7, we get 
a value for the critical removal of qc — 0.988(2). The 
corresponding rate flsphgw is ^sphgw — 2.61(7). 



VI. CONCLUSIONS 

In summary, we have explored how correlations be- 
tween neighboring edge weights can occur in random net- 
works. The high (low) weights can appear concentrated 
in certain areas of the graph, a configuration that has 
a considerable effect on transport properties. To study 
this phenomenon, we have proposed a simple method to 
introduce weight correlations in otherwise uncorrelated 
graphs. These models show that weight correlations can 
appear independently of any other property of the net- 
work, although they could be also coupled to some char- 
acteristic of the vertices such as the degree, hidden vari- 
ables, etc. This method allow us to study, not only quali- 
tatively but also quantitatively, the type and intensity of 
these correlations. Leading us to test several estimators: 
CTu,, the generalized disparity and the range p, being the 
latter the best of the three. 

Once we found a tool to measure the intensity of weight 
correlations, we have focused on how the transport prop- 
erties of the network become affected by these correla- 
tions. The so called superhighways of our model 
have been studied as a function of the intensity of weight 
correlations. The conclusion, that seems to be generaliz- 
able to other networks, is that stronger (positive) correla- 
tions imply bigger and weightier superhighways, improv- 
ing thus the performance of the network on transport in 
orders magnitude. 

Finally, we have also considered data from two real- 
world networks, a collaboration graph and a transporta- 
tion (airports) network. Both cases present positive 
weight- weight correlations. The results on their super- 
highways, ^sphgw, also prove that weight correlations are 
without doubt an important factor to take into account 
in the study of transport on real networks. 
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